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We have computed o? F^s for the hole-doped cuprates within the framework of the one-band Hub- 
bard model, where the full magnetic response of the system is treated properly. The d-wave pairing 
weight a^Fd is found to contain not only a low energy peak due to excitations near (7r,7r) expected 
from neutron scattering data, but to also display substantial spectral weight at higher energies due to 
contributions from other parts of the Brillouin zone as well as pairbreaking ferromagnetic excitations 
at low energies. The resulting solutions of the Eliashberg equations yield transition temperatures 
and gaps comparable to the experimentally observed values, suggesting that magnetic excitations 
of both high and low energies play an important role in providing the pairing glue in the cuprates. 
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I. INTRODUCTION 

Since the superconducting state in the cuprates evolves 
from the doping of a Mott insulator, it is natural to con- 
jecture that the pairing is driven by magnetic fluctuations 
rather than by phonons. Quantum Monte Carlo (QMC) 
calculations provide evidence for d-wave pairing [l|, [3], 
where the pairing bosons reside predominantly in the 
(transverse) spin channel. Recent debate in this con- 
nection has centered on whether or not the magnetic 
resonance peak is strong enough to account for the con- 
densation energy S 0] ■ Although recent estimates seem 
to be affirmative [5j, they do not take into account com- 
peting pairbreaking effects which enter the Eliashberg 
equations [gJ. Also, there are arguments that high-energy 
excitations play a role0, H, @]. Here we report a com- 
putation of a^F's for the hole-doped cuprates based on 
the one-band Hubbard model, where the full magnetic 
response of the system is included, and the Eliashberg 
equations are then solved selfconsistently to obtain the 
superconducting properties over a wide range of dopings 
and temperatures. The resulting transition temperatures 
and pairing gaps are found to be comparable to experi- 
mental values, showing clearly the viability of the mag- 
netic mechanism in the cuprates. We find that excita- 
tions at both high and low energies are important. 

Early calculations of magnetic pairing in the cuprates 
employed parametrized models of the susceptibility. The 
analysis of Radtke, et al. (RULN)[i3l invokes neu- 
tron scattering measurements, while that of Millis, et al. 
(MMP)[lT| is based on NMR data. The model a^F's so 
obtained lead to divergent predictions concerning the fea- 
sibility of magnetic mechanism Our d-wave pairing 
weight a^Fd contains not only a low-energy peak (LEP) 
from near-(7r,7r) scattering, but also an additional high 
energy feature (HEF) extending to ~1.5 eV dominated 
by other regions of the Brillouin zone (BZ), as well as a 
significant pairbreaking contribution at low energies from 
ferromagnetic fluctuations. The HEF, which was missing 
in the RULN and MMP models turns out to be crucially 



important in producing high transition temperatures and 
pairing gaps. The pairbreaking terms begin to dominate 
as the Fermi energy approaches the van Hove singularity 
(VHS) with increasing doping and can lead to the loss of 
superconductivity. 

Our study bears on the recently discovered 'waterfall' 
or high-energy kink (HEK) features observed over 0.3-1.2 
eV range in the angle-resolved photoemission (ARPES) 
spectra of a number of cuprates. The magnetic sus- 
ceptibility underlying our computation of a^F's yields 
self-energies and dispersions consistent with the water- 
faU effects [13, [l3|, suggesting that the boson responsible 
for the waterfall effects is also a key player in generating 
significant pairing weight in a^Fd and high condensation 
energy in the cuprates. 

The calculations are based on a one-band Hubbard 
Hamiltonian, extended to include pairing interaction. 
Specifically, in terms of susceptibility xo and the Hub- 
bard on-site repulsion [/, we use the singlet pairing 
potential (lij 
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and the mass renormalization potential (Eq. [8] of 
Ref. [ii(a)) 
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where p and p' are the electron momenta, which are con- 
strained to lie on the Fermi surface. Here ]4 [14] is the 
potential contributing to the normal [anomalous] part of 
the self energy. These expressions have been found to give 
transition temperatures in good agreement with QMC 
resultsfi']. The resulting coupling constants in various 
pairing channels a are 
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where V — Vs for the even parity channels. The normal- 
ized weighting function ga — 9a{p) / iNo\vp\) , where Vp 
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FIG. 1: (Color online) Eliashberg functions a'^Fd and a^Fz 
for hole doping x = 0.30 obtained in this work (green line) 
are compared with results of Refs. [l^ (red dotted line) and 
[Til (blue dashed line). 



is the Fermi velocity and Nq = (27r)^ / ga(p)'^<fp/ 
The ga are weighting functions of various symmetry 
of which the most important are the lowest harmon- 
ics of s-wave and dx2_y2 symmetry, with gs = 1 and 
gd = cos(j>xa) — cos{pya). We also define the coupling 
constant Az via the s-wave version of Eq. [3]with V = Vz- 
Then the effective BCS coupling is A^, = Aq/(1 + A^). 
The symmetrized Eliashberg functions then are 
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where V" is the imaginary part of the corresponding V . 

In the presence of strong ma gne tic fluctuations 
Migdal's theorem is not obevedflil. Tl7l|. We have de- 
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veloped a relatively simple approximation scheme [13|, 
which can successfully reproduce the pseudogap and wa- 
terfall phenomena in the normal state of the cuprates 
over the full doping range. In the overdoped regime, this 
scheme reduces to calculating the self-energy in GW ap- 
proximation, using a reduced U = 3.2t and dispersion 
renormalized via = 2. 2l| These values of U and Zq 
yield self-energies in reasonable accord with the QMC 
results [23| and explain the recently observed waterfall ef- 
fects in photoemission spectra of the cuprates. 13j,ll4| We 
therefore expect these parameters to be most appropriate 
near x = 0.27, but to gain some understanding of how 
the band structure would affect superconductivity in the 
absence of pseudogap effects, we solved the Eliashberg 
equations over the full doping range x — 0—0.4, assuming 
Zq and U to be doping independent. A more satisfactory 
procedure would be to let U increase in the underdoped 
regime. However, in the presence of a pseudogap, a ten- 
sor system of Eliashberg equations needs to be solved, 



and that is beyond the scope of the present calculation. 
In short, we proceed thus by solving Eqs. 1 and 2 using 
U = 3.2t and xo renormalized by Zq. In particular, we 
neglect the additional modifications of Migdal's theorem 
in the superconducting state. Despite this limitation, our 
results provide a benchmark for the Eliashberg formula- 
tion in that we do not invoke empirical susceptibilities as 
has been the case in much of the existing literature. 

Concerning technical details, we use a tight-binding 
parametrization of the dispersion of Bi2212, with the bi- 
layer splitting neglected. [l9| xo is first computed within 
the RPA scheme throughout the BZ for frequencies up to 
2.88 eV. a^F's and the A's are then computed from Eqs. 
1-4. Fermi surface restricted Ehasliberg equations 
are finally used to selfconsistently obtain the gap A(a;) 
and renormalization Z{ll!) functions, with Z{0) = Z = 



II. PURE D-WAVE SOLUTION 

Figure [U which compares our typical d-wave pairing 
weights a'^Fd and a'^F^ with RULN and MMP models, 
highlights our key finding. Our a^Fd (green line) in (a) 
displays two clear featuresji^l: A low energy peak (LEP) 
around 40 meV and a broad high energy hump-like fea- 
ture (HEF) extending from ~ 0.5 — 1.0 eV (see also Fig. 
2(a) below). The LEP arises mainly from the magnetic 
response near (7r,7r), but the HEF is connected with the 
response from other parts of the BZ, particularly near 
(7r,0) and (7r/2, ■7r/2). Our LEP in (a) is similar to the 
weights assumed by RULN and MMP. This resemblance 
is not surprising since the RULN model[l3| was designed 
to match neutron scattering data near (7r,7r), while the 
NMR data utilized by MMP IjJ is also most sensitive to 
weight in this part of the BZ. It has long been known that 
neutron scattering near (tt, tt) accounts for only about 
l/8th of the integrated s pec tral weight expected from a 
total scattering sum rule [2^. By basing their estimate 
solely on the neutron scattering data near (7r,7r), RULN 
severely underestimated the total d-wave pairing weight. 
The MMP analysis, based on NMR data, appears to have 
captured more of the weight— although still missing the 
HEF and thus underestimating the total weight. Note 
also from Fig. 1(b) that both models strongly underes- 
timate the renormalization weight a^Fz, which opposes 
the tendency for pairing. 

The negative dip in Fig. 1(a) at energies below 20 meV 
deserves comment. This dip refiects pair-breaking mag- 
netic scattering (PBS) near F and was overlooked in the 
phenomenological RULN and MMP models. For simplic- 
ity, we will refer to these fluctuations as being ferromag- 
netic (FM), although this is strictly so only at F. This 
PBS is related to earlier indications of FM instabilities 
near a VHSjiBI. [2^. A similar scenario of competing d- 
wave pairing vs pairbreaking effects has been discussed 
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FIG. 2: (Color online) Eliashberg functions a^Fd and a^Fz 
over the doping range x — 0.0 — 0.4. Lines of various colors 
refer to different dopings (see legend in (a)). Left hand panels 
(a) and (b) give results over an extended frequency range of 
— 1.5 eV, while right hand panels (c) and (d) highlight the 
low energy region of — 250 meV on an expanded energy scale. 



FIG. 3: (Color online) Doping dependence of: (a) X^; (b) A^j; 
and (c) Ad{T — 0) (left scale) compared to Xd (right scale). In 
(a) three different computations of Xz are compared based on 
the full Vz of Eq. 2 (blue solid line), a simplified Vzo ~ U'^Xo 
(red dotted line), and the estimate N{0)U (green dashed line), 
where A^(0) is the density-of-states at the Fermi energy. 



in the context of electron-phonon pairing 

Figure 2 shows how a'^F's evolve with doping. In 
(a), the pairing weight in the high energy feature of 
a^Fd is seen to increase monotonically with increasing 
doping, displaying an approximate isosbestic point at 
uj ~ 0.24 eV. In the low energy region in (c), position 
of the peak in a^F^ shifts to lower energies with increas- 
ing doping, and the negative pair breaking peak grows 
drainatically, consistent with the suggestion of Kopp et 
al. 26]. The nature of a^Fd is seen to change quite sub- 
stantially as the Fermi energy approaches the VHS at 
around x — 0.39. Interestingly, by comparing (c) and 
(d), the low energy peak in a^F^ is seen to follow that in 
a^Fd to lower energies with doping. 

Figure [3] shows the doping dependence of Xz , Ad and 
the low temperature gap Ad{T = 0). Three different esti- 
mates of Xz are compared in (a) for illustrative purposes. 
Values based on using the bare susceptibility, Vzo — U^Xo 
(red dashed line), are seen to be quite similar to the sim- 
ple estimate N{0)U (green dotted hne), where N{0) is the 
density-of-states at the Fermi energy. The full Vz (blue 
line) on the other hand yields a significant enhancement 
of Xz over that obtained from xo, especially near the re- 
gion of the VHS peak, indicating that the system is close 
to a magnetic instability. Note that is positive for 
dopings less than w 0.4, but as the Fermi energy enters 



the region of the VHS with increasing doping. Ad rapidly 
becomes large and negative due to FM fluctuations, (c) 
shows that this doping dependence of Ad is well correlated 
with that of the pairing gap. We stress that these results 
hold for a pure dx^-y^ order parameter. Harmonic con- 
tent plays an important role, as will be discussed below, 
Section III. 

We turn now to discuss our solutions of the Eliashberg 
equations. Following common practice, we proceeded by 
discretizing the a^F^s on the real frequency axis. [2^ We 
find that our results are sensitive to the number of 
points in the mesh. For the present calculations, based on 
a 768-point non- uniform mesh over — 2.88 eV, the gap 
A(a;) is approximately converged in the low-cj regime, al- 
lowing us to extract Ad{T). Figure 4(a) shows typical re- 
sults for the real part of A(cl') for a range of temperatures 
at X = 0.10. The prominent oscillations in A(w) curves 
are the well-known consequence of discretizing a^F's in 
solving the Eliashberg equations. [2^ We define the gap 
by taking the intersection of the A(w) = to line (thin 
black line in Fig. 4(a)) with the A(cj) curve. 

Fig. 4(b) shows how the computed low-energy gap Ad 
evolves with temperature at various dopings. Due to the 
difficulty of finding well-converged solutions when A is 
small, we calculate Ad(r) at a few low temperatures, and 
use a fit to a c?-wave BCS gap to estimate Tc- We find 
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FIG. 4: (Color online) (a) Real part of the gap function 
A{uj,T) at doping a;=0.10 as a function of frequency for a 
series of temperatures T (see legend). Thin black line is the 
plot of A = tj used to obtain the low energy gap as discussed 
in the text, (b) Computed temperature dependence of the 
low-energy gap Ad{T) at various dopings x (see legend). 
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FIG. 5: (Color online) Typical computed SIS tunneling spec- 
trum at a; = 0.10. Inset shows the high energy tail on an 
expanded scale. 



2Ad(0)/fcsTc ~ 3.2 for different dopings. The resulting 
Tc's are somewhat smaller than QMC valuesj^j, perhaps 
the effect of a finite t'. 

It is striking that the gap features in Fig. 4(a) ex- 
tend to very high energies, raising the obvious ques- 
tion as to how this high-energy tail would show up in 
tunneling spectra, [s^ Insight in this regard is provided 
by Fig. 5, where we show a typical tunneling spectrum 



computed SlJ within our model. [Tunneling spectra com- 
puted at other dopings are similar, except that the fea- 
tures scale with A^^.] The weight in Fig. 5 at energies 
above the peak-dip-hump feature is seen to be quite small 
with weak energy dependence (see inset) and would not 
be readily observable in the presence of an experimental 
background. 



III. LOW VS HIGH ENERGY PAIRING GLUE 

Within the present model, the LEP and HEF both play 
an essential role in generating large gaps. For example, 
at a; = 0.3, the HEF by itself produces a gap of only ^ 
0.4 meV, while the LEP is virtually nonsuperconducting, 
even though the full a^F^ yields a gap of 5.5 meV. [To 
be definite, we separate a^F into LEP and HEF at the 
minimum in a^F, ujmin = 0.3 eV.] Similarly, for a;=0.1, 
LEP [HEF] by itself has a ~3 [0.4] meV gap, with a 
combined gap of ~17 meV, with uomin = 0.16 eV. 

This behavior can be readily understood from a 2-A 
model. [13] Since this is a purely electronic mechanism, 
we use a modified Allen-Dynes formula 33|, |34 ) 



Ic- ^^^exp[ )- ^_^exp[ ), (5) 



A(0) = 3MTc, with 



a^F{Lo) 



and 



2 f°° 
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UJ 



(6) 
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The Allen-Dynes equation has a well-known 
limitation 3J] that it predicts a maximum Tc — a-';„/1.2, 
whereas the Eliashberg equations have a solution that 
grows without limit ~ as A ^ oo. We find that 
this leads to an underestimate of A^^p, while the 
model provides good estimates for the remaining gaps. 
For instance, for x = 0.1, Xlep = AijEF=0.15, 
^in,LEP = 83 meV, uJin,HEF = 530 meV, so 
^LEF = 0.26 meV, and /^hef = 1-4 meV. When 
both features are combined, aj;„ = 200 meV and 
Ad = 0.3, leading to Ad = 19 meV, in good agreement 
with the full calculation. While the Allen-Dynes model 
is highly simplified, it does capture the observed trend 
that both peaks contribute significantly. Physically, the 
effective A is in the weak coupling regime, A ~<< 1, 
so high Tc arises from the large u)im and the large 
boost from combining LEP and HEF arises since 
g-i/2A 2e~^^^. Clearly, an electron-phonon coupling 
could play a similar role in further enhancing Tc. 
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FIG. 6: (Color online) Doping dependence of: Xd (blue 
line), As (red line), Xdxy (green line), Ap (violet line), and 
Asp (brown line), calculated from a 15x15 harmonic matrix 
in each symmetry sector. 



IV. COMPETING ORDER PARAMETER 
SYMMETRIES 

The above calculations have been limited to a pure 
d-wave gap symmetry, without harmonic content. In 
tetragonal symmetry there are five symmetry classes of 
superconducting gap, and each class can involve higher 
harmonics of the given symmetry [s^. While we have not 
solved the tensor Eliashberg equations, it is straightfor- 
ward to generalize the A calculations to include harmonic 
structure and to calculate the leading A eigenvalue for 
each symmetry class. The results are shown in Fig. |6] fol- 
lowing the analysis of Ref.[36|. We see that: (1) The pure- 
ed analysis of Section II holds in the low doping regime; 
(2) Near the VHS, harmonic content stabilizes d-wave 
symmetry, leading to the largest gaps; (3) In this regime, 
other symmetries can become comparable to d-wave. In 
particular, there is a tendency toward s-wave pairing in 
the overdoped case. 



The scarcity of high- Tc superconductors arises in part 
from the fact that when superconducting pairing is suf- 
ficiently strong, corresponding and stronger instabilities 
arise in other channels. We have for the most part ne- 
glected the effects of competing phases, but it is clear that 
they will be significant, both in the underdoped regime 
and near the VHS. Near the VHS pairbreaking ferromag- 
netic scattering increases sharply, strongly suppressing a 
pure d-wave gap. While we find that inclusion of har- 
monic content could stabilize a d-wave superconductor 
even at the VHS, we have not accounted for a competing 
FM instability. Indeed, Storey, et al. [3§| find that in Bi- 
2212 the VHS induces strong pair breaking, suppressing 
superconductivity, so that optimum Tc falls at a doping 
below the VHS. This is consistent with the evidence for 
strong FM pairbreaking adduced by Kopp et al.[2^ 

To conclude, we have demonstrated that, when real- 
istic a^F^s are used to solve the Eliashberg equations, 
the magnetic mechanism is capable of producing transi- 
tion temperatures and pairing gaps comparable in size to 
those observed experimentally in the cuprates. The low 
values of these key superconducting properties found in 
earlier calculations are directly attributable to the fact 
that neutron scattering sees only a fraction of the total 
magnetic spectral weight in these materials. 
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V. CONCLUSIONS 

In summary, we have shown that within the present 
model the d-wave pairing weights a^Fd and a^F^ extend 
to very high energies of ~ 1 eV when the magnetic re- 
sponse of the system is properly taken into account. The 
associated superconducting gap is quite substantial, be- 
ing around 16 meV at low dopings. a^Fd is found to 
contain not only the expected low energy peak (LEP) 
below 200 meV, but a previously unrecognized high en- 
ergy feature (HEF) over 0.3-1.2 eV. We find that the LEP 
and HEF both play an important role in yielding a large 
gap in our model. This suggests that electron-phonon 
coupling could be important for further enhancing Tc as 
suggested by the isotope effect [sT]. The gap spectrum 
A(w) generally extends to the limit of a^F (^^1.5 eV), 
and may provide insight into a number of anomalous fea- 
tures connected with optical properties of the cuprates, 
summarized for example in Ref . l38l 
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